Non-unital non-Markovianity of quantum dynamics 
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Trace distance is available to capture the dynamical information of the unital aspect of a quantum 
process. However, it cannot reflect the non-unital part. So, the non-divisibility originated from the 
non-unital aspect cannot be revealed by the corresponding measure based on the trace distance. We 
provide a measure of non-unital non-Markovianity of quantum processes, which is a supplement to 
Breuer-Laine-Piilo (BLP) non-Markovianity measure. A measure on the degree of the non-unitality 
is also provided. 
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I. INTRODUCTION 



Understanding and characterizing general features of 
the dynamics of open quantum systems is of great im- 
portance to physics, chemistry, and biology [1]. The non- 
Markovian character is one of the most central aspects of 
an open quantum process, and attracts increasing atten- 
tions |2l-[l6|. Markovian dynamics of quantum systems 
is described by a quantum dynamical semigroup [1], [l7| , 
and often taken as an approximation of realistic cir- 
cumstances with some very strict assumptions. Mean- 
while, exact master equations, which describe the non- 
Markovian dynamics, are complicated [9]. Based on the 
infinitesimal divisibility in terms of quantum dynamical 
semigroup. Wolf et al. provided a model-independent 
way to study the non-Markovian features (H, Later, 
in the intuitive picture of the backward information flow 
leading to the increasing of distinguishability in interme- 
diate dynamical maps, Breuer, Laine, and Piilo (BLP) 
proposed a measure on the degree of non-Markovian be- 
havior based on the monotonicity of the trace distance 
under quantum channels [4], as shown in Fig. [TJ The 
BLP non-Markovianity has been widely studied, and ap- 
plied in various models [18-23J. 

Unlike for classical stochastic processes, the non- 
Markovian criteria for quantum processes is non-unique, 
and even controversial. First, the non-Markovian crite- 
ria from the infinitesimal divisibility and the backward 
information flow are not equivalent [19, 20]. Second, sev- 
eral other non-Markovianity measures, based on different 
mechanism like the monotonicity of correlations under 
local quantum channels, have been introduced [T3| . 
Third, even in the framework of backward information 
flow, trace distance is not the unique monotone distance 
for the distinguishability between quantum states. Other 
monotone distances on the space of density operators can 
be found in Ref. [24] , and the statistical distance (25|, |26| 
is another widely-used one. Different distance should not 



be expected to give the same non-Markovian criteria. 
The inconsistency among various non-Markovianity re- 
flects different dynamical properties. 

In this paper, we show that the BLP non-Markovianity 
cannot reveal the infinitesimal non-divisibility of quan- 
tum processes caused by the non-unital part of the dy- 
namics. Besides non-Markovianity, "non-unitality" is an- 
other important dynamical property, which is the neces- 
sity for the increasing of the purity Trp^ under quan- 
tum channels [l^l and for the creating of quantum dis- 
cord in two-qubit systems under local quantum chan- 
nels [28]. In the same spirit as BLP non-Markovianity, 
we define a measure on the non-unitality. As BLP non- 
Markovianity is the most widely used measure on non- 
Markovianity, we also provide a measure on the non- 
unital non-Markovianity, which can be conveniently used 
as a supplement to the BLP measure, when the quantum 
process is non-unital. We also give an example to demon- 
strate an extreme case, where the BLP non-Markovianity 
vanishes while the quantum process is not infinitesimal 
divisible. 

This paper is organized as follows. In Sec. Ull we give 
a brief review on the representation of density operators 
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Figure 1: (Color online) Sketch of the information flow picture 
for non-Markovianity [4]. According to this scenario, the loss 
of distinguishability of the system's states indicates the infor- 
mation flow from the system to the reservoir. If the dynamics 
is Markovian, the information flow is always outward, repre- 
sented by the green thick arrow. Non-Markovian behaviors 
occurs when there is inward information flow, represented by 
the orange thin arrow, bringing some distinguishability back 
to the system. 
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and quantum channels with Hermitian orthonormal op- 
erator basis, and various measures on non-Markovianity. 
In Sec. nil) we investigate the non-unitahty and the 
non-unital non-Markovianity and give the corresponding 
quantitative measures respectively. In Sec. IIV| we apply 
the non-unital non-Markovianity measure on a family of 
quantum processes, which are constructed from the gen- 
eralized amplitude damping channels. Section |Vl is the 
conclusion. 



II. REVIEW ON QUANTUM CHANNELS AND 
NON-MARKOVIANITY 

A. Density operators and quantum channels 
represented by Hermitian operator basis. 

The states of a quantum system can be described by 
the density operator p, which is positive semidefinite and 
of trace one. Quantum channels, or quantum operations, 
are completely positive and trace-preserving (CPT) maps 
from density operators to density operators, and can be 
represented by Kraus operators, Choi-Jamiolkowski ma- 
trices, or transfer matrices (29l-[32|. 

In this work, we use the Hermitian operator basis to 
express operators and represent quantum channels. Let 
{A^ I /i = 0, 1, • • • , — 1} be a complete set of Hermitian 
and orthonormal operators on complex space C^, i.e., 
satisfies = A^ and (A^,A^) := Tr(A],Ai.) = S^j,. Any 
operator O on can be express by a column vector 
r := (ro, ri, • • • , r^2_i)T through 

0=Y.r,{0)\, (1) 

/i=0 

with r^{0) := (A^,0). Every r^{0) is real if O is Her- 
mitian. 

In the meantime, any quantum channel 8: p ^ S{p) 
can be represented by T{S) : r{p) r[£{p)] via 

r[£{p)] = T{£)r{p), (2) 

where T{S) is a x (P real matrix with the elements 

V(^):=(A^,f(A,)). (3) 

Furthermore, one can easily check that 

T{£io£2)=T{£i)T{£2) (4) 

for the composition of quantum channels. Here £i o £2 
denotes the composite maps £i{£2{p))- 

Taking into the normalization of the quantum states, 
i.e., Tr(p) = 1, ro can be fixed as ro(p) = 1/Vd for any 
density operator p by choosing Aq = t/^/d with 1 the 
identity operator. In such a case, A^ for /i = 1, 2, • • • ^cP — 
1 are traceless and generate the algebra sn{d). This real 
parametrization r^{p) for density operators is also called 



as coherent vector, or generalized Bloch vector (33-[36j. 
In order to eliminate the degree of freedom for the fixed 
ro, we use the decomposition r = (ro,r)^. Therefore, 
any density operator p can be expressed as 

P=-,+r-X, (5) 

with r the generalized Bloch vector and A represents 
(Ai, A2, • • • , A(^2_x)^. Under this frame, quantum chan- 
nels can be represented by the affine map 113 

r{£{p)) = M{£)r{p)+c{£), (6) 

where M{£) is a real matrix with the dimension d? — 1 
and the elements of the vector c{£) reads 

[c(f)]^ = (A^,^(l))M (7) 

for /i = 1, 2, • • • , - 1. Comparing Eq. ^ with Eq. (O, 
one could find that 

T^,{£) = [M{£)U, (8) 

for /i, = 1, 2, • • • — Thus, T{£) can be decomposed 
into the following sub-blocks: 



" 1 


Olx(d2-l) 


\/~dc 


M 



Reminding that a quantum channel £ is said to be 
unital if and only if ^^(Jl/d) = l/d [37], one could find 
that the necessary and sufficient condition for a unital 
map is that c{£) = 0, namely, 

c{£) = <^ £ is unital. (10) 

Thus, c{£) describes the non-unital property of the quan- 
tum channel £. The necessary and sufficient condition 
above could be easily proved by realizing that the Bloch 
vector of U/d is zero vector, i.e., r = 0. Based on the 
sub-block form of T(^), c{£) = is equivalent to that 
T{£) is block diagonal, i.e., T{£) = diag(l, M(£^)). 

Whether a quantum channel £ is completely positive 
(CP) can be reflected by the Choi-Jamiolkowski ma- 
trix [13, m 

C{£):={£^t){\n){n\), (11) 

where \Q) = ^ Xlj=o b)^b) maximally entangled 

state. Here {\j)} is a basis in Hilbert space. £ is CP if and 
only if the Choi-Jamiolkowski matrix is positive. With 
the Hermitian operator basis, is di (P x dP matrix 

and can be written in the form [33l | 

Substituting this formula into Eq. (pTj) and utilizing 
Eq. ([3]), one could express the Choi-Jamiolkowski matrix 
as 

(f-i 

C{£) = ^T^^i^)^^ ® (13) 
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If £ is unital, it can be reduced into 

C(f) = ^ 1 1,2^,2 +d I . (14) 

B. Non-divisibility and non-Markovianity 

Without the presence of correlation between the open 
system and its environment in the initial states, the re- 
duced dynamics for the open system from t = to any 
t > can be expressed as 



St^Q-p^TlE [U(t) {p®pE)U(t)^] 



(15) 



which is a quantum channel. This indicates that Et^Q 
is CPT. The unitary operator U{t) describes the time 
evolution of the closed entirety, and p^; is the initial 
state of the environment. A quantum process £t •= ^t,o 
is said to be infinitesimal divisible, also called as time- 
inhomogeneous or time-dependent Markovian, if it satis- 
fies the following composition law [2] 



(16) 



for any t2 >ti> 0, where £t2,ti is also completely posi- 
tive and trace preserving. 

Various measures on the degree of the non-Markovian 
behavior of quantum processes have been proposed and 
investigated [3, flll - [l3| . Almost all of the measures on 
the non-Markovianity can be classified into three kinds, 
base on the degree of the violation of the following prop- 
erties owned by the infinitesimal divisible quantum pro- 
cess: 

(i) Monotonicity of distance D under CPT maps. That 
is D{£{pi),£{p2)) < D{pi^p2) for any quantum chan- 
nel f, where D{pi^p2) is an appropriate monotone dis- 
tance under CPT maps on the space of density opera- 
tors [24], including trace distance, Bures distance, sta- 
tistical distance, relative entropy, and fidelity (although 
fidelity itself is not a distance, it can be used to con- 
struct monotone distances) and so on. Some measures 
on non-Markovianity by increasing of the monotone dis- 
tance during the mediate dynamical maps ft2,ti ^^^^ 
been given and discussed in Refs. (3, [l2j. 

The typical measure of this type, which would be used 
later in this paper, was first proposed by Breuer, Laine, 
and Piilo in Ref. ^], based on the monotonicity of trace 
distance 1361, l37l 



Ar(/0i,P2) := 2Tr|pi 



P2\ 



(17) 



where \0\ := V 0^0. Interpreting the increase of the 
trace distance during the time evolution as the informa- 
tion flows from the environment back to the system, the 
definition of the BLP non-Markovianity is defined by 



ArBLp(^t) := max / dt a (t, pi, /)2) 



where 

a(t,pi,p2) := ^Ar(pi(t),P2(t)), (19) 

and pi{t) = ft(pi) for z = 1, 2 are two evolving states. 

(ii) Positivity of the Choi-Jamiolkowski matrix for 
CPT maps. The Choi-Jamiolkowski matrix C{£) > 
if and only if f is a quantum channel, namely, f is a 
CPT map. Some measures on non-Markovianity by the 
negativity of the Choi-Jamiolkowski matrix for mediate 
dynamical maps 5^2, ti^^^^ been given and discussed in 
Refs. 

In this work we would use one of these measures, 
which was proposed by Rivas, Huelga and Plenio (RHP) 
in Ref. [6]. They utilize the negativity of the Choi- 
Jamiolkowski matrix C for the mediate dynamical maps 
with the definition 



MiHp(^t) := / g{t)dt, 

Jo 



(20) 



where 



(18) 



,it) := lim - \ (21) 

(iii) Monotonicity of correlations E under local quan- 
tum channels. That is E [{S (g) 1) (p^^)] < E{p^^) for 
any local quantum channel 8, where E is an appropriate 
measure for the correlations in the bipartite states p^^, 
including entanglement entropy and the mutual informa- 
tion. The corresponding measures on non-Markovianity 
are given and discussed in Refs. j6L Il3|. 



III. NON-UNITAL NON-MARKOVIANITY 

The non-Markovianity measure A/rlp is available to 
capture the non-Markovian behavior of the unital aspect 
of the dynamics. But for the non-unital aspect, it is not 
capable. To show this, we use the Hermitian orthonormal 
operator basis to express states and quantum channels. 
Utilizing Eq. ([5]), the trace distance between two states 
pi and p2 is given by 



Dtr{pi,P2) = -Tr|[r(pi)-r(/,2)]-A|. (22) 
Therefore, for the two evolving states, we get 



D^Api{t),P2{t)) = -Tr|M(£:t)[r(pi) -r(p2)] • A|, (23) 

where pi, p2 are initial states of the system. 

From this equation one can see that the trace distance 
between any two evolved states is irrelevant to the non- 
unital part c{St) of the time evolution. Then if there 
are two quantum channels, whose affine maps are r 
Mr + ci and r i-^ Mr + C2, respectively, the characteristic 
of trace distance between the evolving states from any 
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two initial states cannot distinguish these two channels. 
More importantly, c{£t) may cause the non-divisibility of 
the quantum process f^, and this cannot be revealed by 

A^BLP. 

On the other hand, the non-unital part c{St) has its 
own physical meaning: c{St) ^ is necessary for the 
increasing of the purity V{p) = Tr(p^) [27]. In other 
words. 



c{St) = 



(24) 



Besides the non-Markovian feature, the non-unitality is 
another kind of general feature of quantum processes. In 
analogy to the definition of BLP non-Markovianity, we 
defined the following measure on the degree of the non- 
unitality of a quantum process: 



NnuiSt 



max 

Po 



1 




/^T'(£,(po))>0 


At 



dt, (25) 



where po is the initial state. Obviously, Nnu{£t) vanishes 
if c(ft) = 0. 

Since the non-unital aspect of the dynamics, which is 
not revealed by the trace distance, has its own speciality, 
we aim to measure the effect of non-unitality on non- 
Markovian behavior. However, a perfect separation of 
the non-unital aspect from the total non-Markovianity 
may be infeasible. Therefore we require a weak version 
Afnu for measuring non-unital non-Markovianity to sat- 
isfy the following three conditions: (i) A^nu vanishes if 
£t is infinitesimal divisible, (ii) Afnu vanishes if £t is uni- 
tal, (iii) A/'nu should be relevant to c{£t). Based on these 
conditions, we introduce the following measure 



max 



nu>0 



(26) 



where X := {gr | < r < oo} with := Sr{'^/d) is the 
set of the trajectory states which evolve from the maxi- 
mally mixed state, and 

(Tnnit^gr) := [ft(^o), (^r)] , (27) 

with D{pi,p2) an appropriate distance which will be dis- 
cussed below. The first condition is guaranteed if we 
require that D is monotone under any CPT maps, i.e., 
D[8{pi),£{p2)] < D{pi,p2) for any quantum channel £. 
For the unital time evolution, the set X = {t/d} only 
contains the maximally mixed state, so the above de- 
fined A/'nu vanishes, and the second condition is satisfied. 
The third condition excludes the trace distance. 

In this paper, we use the Bures distance which is de- 
fined as 



1^b(pi,P2) = V2 [l-F{p 



(28) 



where 



F{pi,P2) = Tv\^^\ = T:tJ^P2Vp^ (29) 



is the Uhlmann fidelity 
10 



between pi and p2. Here 



V 0^0. Bures distance is an appropriate distance 
for A/'nu because it obeys the monotonicity under CPT 
maps [24] and is relevant to c{£t). As here only the 
monotonicity of distance is relevant, for simplicity, we 
can also take the square of the Bures distance or just the 
opposite value of Unlmann fidelity as a simple version of 
monotone "distance" [12]. Quantum relative entropy 
S{pi\\p2) = Tr [pi(lnpi — lnp2)], or its symmetric ver- 
sion 6'sym(pi IIP2) := S{pi\\p2)^S{p2\\pi), is another qual- 
ified candidate for the distance. Noting that when the 
support of pi is not within the support of p2, namely, 
supp(pi) ^ supp(p2), S{pi\\p2) will be infinite, so in such 
cases, quantum relative entropy will bring singularity to 
the measure of non-Markovianity. Also, Hellinger dis- 
tance [41] is qualified. Although all of these distances 
are monotone under CPT maps, they may have different 
characteristics in the same dynamics, see Ref. [42]. 

The difference between non-unital non-Markovian 
measure defined by Eq. (|26]) and the BLP-type measures, 
including those which use other alternative distances, is 
the restriction on the pairs of initial states. Comparing 
with the BLP-type measures relying on any pair of initial 
states, the non-unital non-Markovianity measure only re- 
lies on the pairs consisting of the maximally mixed state 
and its trajectory states. On one hand, this restriction 
makes the non-unital non-Markovianity measure vanish 
when the quantum processes are unital, no matter they 
are Markovian or non-Markovian; on the other hand, 
this restriction refiects that non-unital non-Markovianity 
measure reveals only a part of information concerning the 
non-Markovian behaviors. 



IV. EXAMPLE 

To illustrate the non-unital non-Markovian behavior, 
we give an example in this section. We use the general- 
ized amplitude damping channel (GADC) as a prototype 
to construct a quantum process. The GADC can be de- 



scribed by £{p) 
{Ei} given by [37 



EipE\ with the Kraus operators 



^1 = VP 



E2 



E^ 




(30) 

(31) 
(32) 
(33) 



where p and 77 are real parameters. Note that for any 
p G [0,1] and any r] G [0,1], the corresponding ^ is a 
quantum channel. For a two-level system, the Hermitian 
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orthonormal operator basis can be chosen as A = cr 
where cr — {(Jx^cTy^az} is the vector of Pauh matrices. 
With the decomposition in Eq. (|5]) , the affine map for the 
Bloch vector is given by r{S{p)) ^ M(f )r(p) + c(f ) 
where 



(a) 



M{S) 




(34) 



(35) 



The GADC is unital if and only if p = 1/2 or = 1. 
When 7^ = 1, M(^) = 1, the map is identity. 

A quantum process can be constructed by making the 
parameter p and r] to be dependent on time t. For sim- 
phcity, we take pt = cos^ ujt and rjt = e~^, where uj is 
a constant real number. This is a legitimate quantum 
process, because St is a quantum channel for every t > 0, 
and £t=o is the identity map. 

First, let us consider the A/blp for this quantum pro- 
cess. For any two initial states pi and p2, we have the 
trace distance 



MiStMp,)-rip2)] 



V2 



V2 



\M{£tMp,)-v{p2)]\, (36) 



where |r| = y'r • r is the Euclidean length of the vector 
r, and we used the equality 



(a • a){h • tr) = (a • b)lL + ztr • (a x b) 



(37) 



for Pauli matrices. Denoting t{pi) — r(/>2) by (x,y,z)^, 
we get 



Ar [St{pi),St{p2)] 



-t/2 



V2 



(38) 



which implies -^Dt^ [St{pi), £t{p2)] ^ for every time 
point t > and for any real numbers and z. Thus, 
the BLP non-Markovianity vanishes, i.e., A/blp (^t) = 0, 
although £t may be not infinitesimal divisible, which will 
become clear later. 

In order to investigate whether £t is infinitesimal divis- 
ible or not, we shall apply A/'nu in the above model. The 
trajectory of the maximally mixed state under St reads 



St{,o) = -t^Cf^ = - 



l^Wt 
l-Wt 



(39) 



where 
Wt :-- 



(2pt-l)(l-r?t)=cos(2wt)(l 



(40) 



Bures distance 
Trace distance 




0.5 



1.0 

t 



1.5 



2.0 




Figure 2: (Color online) (a) Evolution of trace distance and 
Bures distance between two evolving states of a two-level sys- 
tem under the variant generalized amplitude damping chan- 
nel, initially from the maximal mixed states = 1/2 and its 
trajectory state Qt = ^t(^o), respectively, (b) The evolution 
of g{t) defined by Eq. ()2ip , whose integral with respect to time 
t is RHP measure for non-Markovianity. In these figures, the 
parameters are taken as r = 10 and = 5. 



the corresponding evolving states: 
StiQr) = ^1L + (M,C, + Ct)~ 

^ 1 i^Wt^mWr 

2 I l-Wt- VtWr 




Then the fidelity reads 



where 



F[St{go),St{gr)] = -{h^^h^), 



h+ := ^{l^Wt){l^Wt^r]tWr), 
h- := V{^-Wt){l-Wt-VtWr). 



(43) 

(44) 
(45) 



To compare with the behavior of trace distance, we also 
get Ar [St{go),St{gr)] = IVtWrl /2. With the expressions 
r]t = and pt = cos^ cjt, it is 



Taking these trajectory states as the initial states, we get 



Ar [St{go),St{gT)] 



|cos2a;r| (1 - e"^). (46) 
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In Fig. [2Ja), we can see that while the trace dis- 
tance between the evolving states St{go) and StiQr) 
monotonously decreases with the time the Bures dis- 
tance increases during some intermediate time intervals. 
From Eq. (|46|) . one can see although Dtr[St{go),St{gT)] 
depends on Wr , it does not depend on Wt . Actually, from 
Eq. ([38|) one could find that for any two initial states, the 
trace distance between the evolving states is independent 
on Wt. In this sense, the BLP non-Markovianity treats 
a family of quantum processes, which only differ with p^, 
as the same one. Meanwhile, A/'nu reveals the effects of 
Pt on the infinitesimal non-divisibility and is capable of 
measuring it. 

In order to compare with BHP measure, we also cal- 
culate the g{t) defined by Eq. (|2l]). We get 

g{t) = l[\i-f{t)\^\f{t)\-i] (47) 

with 

f{t) := -ujsm{2ujt) (l - e"^) +cos^(cjt). (48) 

The mediate dynamical maps St-\-e,t with infinitesimal 
e are not completely positive when g{t) > 0. From 
Fig. El^b), we can see that the increasing of the Bures 
distance occurs in the regimes where g{t) > 0, which co- 
incides with the monotonicity of Bures distance under 
CPT maps. 

V. CONCLUSION 

In conclusion, we have shown that the measure for 
non-Markovianity based on trace distance cannot reveal 



the infinitesimal non-divisibility caused by the non-unital 
part of the dynamics. In order to refiect effects of the 
non-unitality, we have constructed a measure on the 
non-unital non-Markovianity, and also defined a mea- 
sure on the non-unitality, in the same spirit as BLP non- 
Markovianity measure. 

Like non-Markovianity, the non-unitality is another in- 
teresting feature of the quantum dynamics. With the de- 
velopment of quantum technologies, we need novel the- 
oretical approaches for open quantum systems. It is ex- 
pected that some quantum information methods would 
help us to understand some generic features of quantum 
dynamics. We hope this work may draw attention to 
study more dynamical properties from the informational 
perspective. 
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